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This Letter reports a measurement of the cross section for producing pairs of central prompt 
isolated photons in proton-antiproton collisions at a total energy ^/s = 1.96 TeV using data cor- 
responding to 9.5 fb~^ integrated luminosity collected with the CDF II detector at the Fermilab 
Tevatron. The measured differential cross section is compared to three calculations derived from the 
theory of strong interactions. These include a prediction based on a leading order matrix element 
calculation merged with parton shower, a next-to-leading order, and a next-to-next-to-leading order 
calculation. The first and last calculations reproduce most aspects of the data, thus showing the 
importance of higher-order contributions for understanding the theory of strong interaction and 
improving measurements of the Higgs boson and searches for new phenomena in diphoton final 
states. 
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lisions is a significant, irreducible background in searches 
for a low-mass Higgs boson decaying into a photon pair 
as well as in searches for new phenomena, such as 
extra spatial dimensions and two-body Q or cas- 

cade [6| decays of new heavy particles. Precise measure- 
ments of the production cross sections for diphotons as 
functions of various kinematic variables and their theo- 
retical understanding are important for these searches. 
The better the prompt diphoton background is under- 
stood, the smaller uncertainties are introduced in these 
searches. After the recent discovery of the Higgs boson- 
like particle at the LHC , a better understanding of the 
background is important for improvements in the preci- 
sion of the measurements of the production cross sec- 
tion and the decay branching ratio of this particle into 
a photon pair. A precise measurement of the branch- 
ing ratio is of special importance, as this decay pro- 
ceeds through a fermion loop and thus it indirectly con- 
strains the couplings of the Higgs boson-Uke particle to 
fermions, which are more difhcult to extract from di- 
rect decays into fermion pairs. Diphoton production is 
also used to test quantum chromodynamics (QCD), the 
theory of strong interaction, both in the perturbative 
scheme (pQCD), which is a good approximation at high 
energies, and in non-perturbative schemes, such as soft- 
gluon resummation methods, which provide important 
corrections in certain lower-energy kinematic regions [8j . 
Diphotons are expected to be dominantly produced by 
quark-antiquark annihilation qq 77 and, in kinematic 
regions where gluons dominate the parton distribution 
functions (PDF) , by gluon-gluon fusion gg — > 77 through 
a quark loop amplitude. Prompt photons may also result 
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from quark fragmentation in hard scattering, although a 
strict photon isolation requirement significantly reduces 
the fragmentation contributions. 

Diphoton measurements have been made previously at 
fixed-target and collider experiments (lol- 12|. Re- 



cent measurements have been made both at the Teva- 
tron [H, 14| and at the LHC 15 1, which offer a consis- 
tent picture on the accuracy and limitations of the the- 
oretical calculations in reproducing the data. The AT- 
LAS measurement [l^ found diphoton production fea- 
tures in proton-proton collisions at \/s = 7 TeV analo- 
gous to those observed in proton-antiproton collisions at 
= 1.96 TeV [H,!!!. The most recent CDF measure- 
ment [l3l, using approximately half the full CDF data 
sample, compared the data with pQCD calculations at 
leading order (LO) and next-to-leading order (NLO) in 
the expansion parameter Us, the strong interaction cou- 
pling. Large discrepancies were found between the data 
and a LO matrix-element calculation supplemented with 
a parton shower (PS) model. The inclusion of photons 
radiated from initial- and final-state quarks allowed by 
the shower model substantially improved the agreement 
of the PS calculation with the data. The calculation that 
includes radiated photons was recently used to predict 
the non-resonant background in the search for a low-mass 
Higgs boson decaying into a photon pair using the full 
CDF data set [H. 

This work presents the final diphoton measurements 
from CDF using the fuU data set collected in 2001-2011 
corresponding to a total integrated luminosity of 9.5 fb~^. 
The results are compared with all the available state-of- 
the-art calculations under a variety of kinematic condi- 
tions 17|, including an improved set of calculations not 



discussed in the previous work [14 1. 

The reported measurement is using data collected with 
the Collider Detector at Fermilab (CDF) [ISl, at the 
Tevatron pp collider. The CDF detector includes a cen- 
tral spectrometer inside a 1.4 T axial magnetic field, sur- 
rounded by electromagnetic and hadronic calorimeters 
and muon detection chambers. The inner spectrome- 
ter measures charged particle trajectories (tracks) with 
a momentum component transverse to the beam (px) 
with a precision of ap^/p^ — 0.07%(GeV/c)^^ The 
pointing-tower-geometry central calorimeters cover the 
region \ri\ < 1.1, with an electromagnetic (hadronic) en- 
ergy resolution of a{ET)/ET = 13.5%/v/lRGeV) © 
1.5% ia{ET)/ET = 50%/V£^T(GeV)® 3%) and a tower 
segmentation of A77 x A(f> ~ 0.1 x 15°, where Et = E sin 9 
is the transverse energy, rj — ~ ln[tan(6'/2)] is the pseudo- 
rapidity, 6 is the polar angle and (j) the azimuth of the 
tower's axis in the coordinate system of the laboratory, 
with polar axis along the proton beam direction and ori- 
gin at the center of the detector. Photons are recon- 
structed in clusters of up to three towers [19| in the cen- 
tral calorimeter only. The pseudorapidity of each photon 
in the event is restricted to the region < 1. A finely- 



segmented detector located at a depth corresponding to 
the maximum development of a typical electromagnetic 
shower measures the energy deposit profile, which is re- 
quired to be consistent with originating from a single 
photon. The photon transverse energy is required to ex- 
ceed 17 GeV for the first photon in the event and 15 GeV 
for the second photon. The transverse energy measured 
by the calorimeter in an isolation cone with radius in 
T] — 4> space of 0.4 around each photon [i^] is required 
not to exceed 2 GeV. 

This measurement employs the same techniques as the 
previous work 1^ . Inclusive diphoton events are selected 
online by requiring two isolated electromagnetic clusters 
with Et > 12 GeV each or two electromagnetic clus- 
ters with Et > 18 GeV and no isolation requirement. 
In the offiine analysis additional requirements are im- 
posed to identify a sample rich in prompt photons. The 
background from events where one or both reconstructed 
photons are misidentified jets is subtracted with a 4 x 4 
matrix technique using the track isolation as the discrim- 
inant between signal and background 2l|. The matrix 
is constructed for each event from the i?T-dependent ef- 
ficiencies of signal and background photons passing the 
track isolation criterion. This technique takes into ac- 
count the full correlations between the two photons in 
the event. An optimal track- isolation threshold of 1 
GeV/c is determined by maximizing the discrimination 
between signal and background Monte Carlo (MC) sim- 
ulation samples. The efficiencies used in this method are 
determined from 7-l-jet and dijet samples generated with 
PYTHIA [22[, subjected to the full detector and online 
event selection simulation [2^, and reconstructed as the 
experimental data. The signal fraction is determined by 
summing the probabilities of an event to be pure signal, 
pure background, and a mixed photon pair. These prob- 
abilities are obtained for each event by multiplying the 
4x4 matrix with the four-dimensional column vector 
of the observation values (0 or 1) for all four combina- 
tions of the first and second photon having track isolation 
larger or smaller than 1 GeV/c and inverting the result. 
The systematic uncertainty in the signal fraction is in the 
range of 15-20%. 

The differential cross section for diphoton production 
is obtained from the histogram of the estimated signal 
yield as a function of each relevant kinematic variable. 
The average cross section in a bin is determined by divid- 
ing the yield by the product of the trigger efficiency, the 
selection efficiency and acceptance, the integrated lumi- 
nosity, and the bin size. The diphoton trigger efficiency 
is derived from data Q. It is consistent with 100% over 
all of the kinematic range with a flat uncertainty of 3%. 
The diphoton selection efficiency accounts for the effects 



from the underlying event from collision remnants 14 1 
and from additional (pile-up) collisions overlapping with 
the collision that produced the photons. The systematic 
uncertainty in the selection efficiency related to the pile- 
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up effect grows linearly from 1.8% for Et < 40 GeV to 
3% for Et = 80 GeV and remains constant above this 
point. A flat 3% uncertainty per photon accounts for 
possible inaccuracies in the PYTHIA model for the un- 
derlying event. This is summed linearly to 6% for two 
photons, since the underlying event is not related with 
prompt photon production and affects only the isolation 
in a way, on the average, symmetric for the two photons. 
A 6% constant uncertainty comes from the integrated lu- 
minosity [i^l . The difference in the photon identification 
efficiency between data and MC is estimated from the 
— >• e+e~ sample [H and added as a systematic uncer- 
tainty to the measurement. The electromagnetic energy 
scale is determined from the mass of the Z° — > e+e~ sig- 
nal. The associated systematic uncertainty is estimated 
to grow linearly from at Et < 40 GeV up to 1.5% at 
Et = 80 GeV and remain constant above this point. All 
systematic uncertainties are added in quadrature. 

In the previous measurement 1J|, the experimental 



results were compared with three theoretical calcula- 
tions: (i) the fixed NLO predictions of the DIPHOX pro- 
gram [25i], including non-perturbative parton fragmen- 
tation into photons at NLO [l^, (ii) the predictions of 
the RESBOS program Q where the cross section is ac- 
curate to NLO, but also has an analytical initial-state 
soft-gluon resummation, and (iii) the predictions of the 
PYTHIA PS program [22| including photons radiated from 
initial- and final-state quarks [IJ]. Within their known 
limitations, all three calculations reproduced the main 
features of the data, but none of them described all as- 
pects of the data. In this Letter, the measurement is com- 
pared with three different calculations: (a) the fixed NLO 
predictions of the MCFM program [27| . including non- 
perturbative parton fragmentation into photons at LO 
[28|, (b) the fixed next-to- next-to-leading order (NNLO) 
predictions of a recent calculation [2^ , and (c) the predic- 
tions of the SHERPA program [s^] , based on a matrix ele- 
ment calculation merged with parton shower (ME-I-PS). 
This calculation features a realistic representation of the 
physics events including initial- and final-state radiation. 
The prediction of MCFM is an alternative calculation to 
DIPHOX, but it has not been tested against any previous 
measurement. The NNLO and SHERPA predictions are 
recent calculations that are expected to reproduce the 
data features better than the previous calculations. 

While the NLO and NNLO matrix elements for dipho- 
ton production include all real and virtual processes at 
fixed order in as, the SHERPA matrix element includes 
only real processes at NNLO. However, by merging the 
matrix element contribution (the hard scattering process) 
with those from the parton shower (cascade radiation 
subprocesses from the initial- and final-state quarks and 
gluons), this calculation accounts for real processes effec- 
tively at all orders in as- It also accounts for some vir- 
tual effects via corrections applied in the parton shower 
subprocesses. The SHERPA calculation is an extension of 



the PYTHIA calculation including photons radiated from 
initial- and final-state quarks which was introduced in 
the previous measurement [l3|- In the default SHERPA 
calculation the scale is adjusted to the event kinemat- 
ics automatically by the program [s^l- An uncertainty 
of this calculation is estimated by the difference from an 
alternative calculation which uses a fixed scale. All cal- 
culations are subject to the experimental kinematic and 
isolation requirements [l7| . Theoretical uncertainties are 
best estimated for the fixed-order NLO and NNLO cal- 
culations, where the scale uncertainties are well-defined. 
The estimation is done by increasing and decreasing the 
scale of each calculation by a factor of two relative to 
the default scale and, for the NLO PDF uncertainties, 
by using the 20 CTEQ6M eigenvectors 

The measured cross section for diphoton production 
integrated over the acceptance is 12.3 ± 0.2stat ± 3.5syst 
pb. The predictions for the integrated cross section are 
12.4±4.4pb from SHERPA, 11.5±0.3 pb from MCFM, and 
11-8^0 g pb from the NNLO calculation. All predictions 
are consistent with the measurement. Figure [T] shows the 
comparisons between the observed and predicted dipho- 
ton distributions in mass M, transverse momentum Pt 
and azimuthal separation A4> between the momenta of 
the two photons in the event. 

All predictions for the mass distribution show a reason- 
able agreement with the data for all calculations above 
the maximum at 30 GeV/c^, particularly in the region 
around M=125 GeV/c^ relevant to measurements of 
the Higgs boson 0. All predictions underestimate the 
data rate around and below the maximum, although the 
NNLO prediction reproduces better the data than the 
other two predictions. The SHERPA prediction tends to 
underestimate the data for M>250 GeV/c^. 

In the Pt spectrum, the MCFM prediction underesti- 
mates the data in the region between 30 and 60 GeV/c, a 
feature also observed in the earlier measurements [13, Hal- 
The other two predictions describe the data fairly well in 
this region. For Pt<20 GeV/c, where soft gluon radia- 
tion becomes important, only the SHERPA prediction pro- 
vides a good description of the data because the parton 
showering provides an effective resummation of multiple 
soft-gluon emission amplitudes. The fixed-order predic- 
tions diverge in the limit of vanishing Pt- The NNLO 
prediction tends to overestimate the data rate for Pt>60 
GeV/c. 

Of special importance is the Ac/) spectrum where all 
PS and NLO predictions examined in the previous pa- 
pers failed to describe the data over the full range. The 
SHERPA model shows the best agreement at larger Acf), 
where the diphoton system acquires substantial trans- 
verse momentum due to multiple soft-gluon emission. 
However, SHERPA progressively underestimates the data 
rate below 1.5 rad. The NNLO calculation is the only 
prediction consistent with the data in the low A(j) tail, 
which contains photon pairs with very low mass and rcl- 
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atively high Pt- This calculation tends to underestimate 
the data rate above 1 rad. The SHERPA and NNLO pre- 
dictions generally are in better agreement with the data 
than MCFM. This shows that higher than NLO contribu- 
tions, included in both calculations in different ways, are 
needed in order to better describe the data. More chan- 
nels open at higher order, such as diphoton production 
associated with the emission of two final-state partons 
(2—^4 channels), which enhance the event rate at high 
Pt and low Ac/). 

The observed cross section enhancements at very low 
diphoton mass (M<30 GeV/c^), moderate diphoton 
transverse momentum (30<Pt<60 GeV/c) and low A(p 
(<1 rad) are correlated. The events involved in this cor- 
relation have a topology of same-side diphotons recoiling 
against at least one hard jet. For some of the contribu- 
tions the cross section is enhanced, such as when the two 
photons are emitted by the same parton and are, there- 
fore, predominantly almost collinear. Enhanced contri- 
butions begin to appear in 2— >-3 subprocesses. The im- 
portance of 2— >3 subprocesses was shown in the previous 
CDF measurement [ij], where the inclusion of photons 
radiated in hard 7-l-jet events substantially improved the 
agreement of the PS calculation with the data with re- 
spect to the simple 2— >-2 diphoton calculation. These 
subprocesses are treated in different ways at different or- 
ders of approximation. At NLO, diphotons emitted from 
the same parton can only appear in the fragmentation 
At NNLO such contributions can re- 



25 



components 

suit directly either from 2— >3 subprocesses, where a quark 
loop is included in the diphoton production amplitude, 
or from tree- level 2— >-4 subprocesses [2^. The SHERPA 
calculation also includes 2— )-4 subprocesses [s^]. Thus 
NNLO and SHERPA describe the observed enhancement 
better than mcfm, which does not include such subpro- 
cesses. 

In summary, the diphoton production cross section, 
differential in kinematic variables sensitive to the parton- 
level processes that govern the reaction, is measured us- 
ing all data collected with the CDF II detector, corre- 
sponding to an integrated luminosity of 9.5 fb~^. This 
measurement is consistent with the past CDF measure- 
ments 12, 3] and supersedes them. The measurement 
uses photons with |?7| < 1 and has sufficiently high preci- 
sion to resolve differences between state-of-the-art theo- 
retical predictions. The results are compared with three 
calculations, which apply complementary techniques in 
predicting the cross section. The NNLO calculation is 
generally consistent with the data, although events with 
very low diphoton mass and high diphoton transverse 
momentum are not accurately described. The ME-I-PS 
SHERPA calculation is also consistent with the data ex- 
cept in the tails of the mass and the low A(j) distributions. 
Both NNLO and SHERPA describe the data better than 
the NLO MCFM calculation in regions sensitive to dipho- 
ton production channels resulting to nearly collinear pho- 



tons. The comparisons show that parton-levcl processes 
of order higher than NLO, which was the standard ap- 
proximation in older calculations, play an important role 
in diphoton production at the current level of experimen- 
tal precision. This conclusion is supported by the findings 
of the recent ATLAS measurement at higher collision en- 
ergy [isl . The inclusion of such processes in background 
calculations is thus important for high precision measure- 
ments of the recently discovered Higgs boson-like particle 
and searches for new phenomena in diphoton final states. 
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Figure 1: The measured differential cross sections for M{-fj) compared with six theoretical predictions discussed 
in the text. The left windows show the absolute comparisons and the right windows show the fractional deviations 
of the data from the theoretical predictions. Note that the vertical axis scales differ between fractional deviation 
plots. 
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Figure 7: The measured differential cross sections for A$(77) compared with six theoretical predictions discussed 
in the text. The left windows show the absolute comparisons and the right windows show the fractional deviations 
of the data from the theoretical predictions. Note that the vertical axis scales differ between fractional deviation 
plots. 
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retical predictions discussed in the text. The left windows show the absolute comparisons and the right windows 
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9 



m „ " CDF II Diphoton 9.5 fb 
-^10 ^ET>15,17GeV, hl<1-0, 
E A R>0.4, lso<2 GeV 

■o 



\ I \ \ ^ \ I \ \ p 

- Data 

DIPHOX CTEQ6M 

.. RESB0SCTEQ6M 

PYTHIAyy+Yi 
.. PYTHIAyy 




a. 



2 2.5 3 
A(|) for P^<M (rad) 



X 

a. 
a 

>< 
a. 
a 



m 
(/) 
cc 

m 

DC 




0.5 1 1.5 2 2.5 3 
A(|) for Pj<M (rad) 



CDF II Diphoton 9.5 fb 

^ET>15,17GeV, hl<1-0, 
A R>0.4, lso<2 GeV 




1.5 2 2.5 3 

A(|) for P.^<M (rad) 



< 

0. 

cc 

CO 

< 

0. 

cc 

CO 



6f 
5! 

4 

3^ 
2 
1 


■d 

5 
4 

3^= 
2 
1 

■I 
5 
4 
3 
2 
1 

-1 



7 



■ Data 
PDF uncertainty 
Scale uncertainty 



3 



0.5 1 1.5 2 2.5 3 

A(|) for Pt<M (rad) 



Figure 9: The measured differential cross sections for A$(77), when Pri'Jj) < -^(77), compared with six theo- 
retical predictions discussed in the text. The left windows show the absolute comparisons and the right windows 
show the fractional deviations of the data from the theoretical predictions. Note that the vertical axis scales differ 
between fractional deviation plots. 



10 



F(77) 



2- 20, 

Q. - CDF II Diphoton 9.5 fb"' 

18-ET>15,17GeV, hl<1-0, 
^ - A R>0.4, IS0<2 GeV 

1 16- 
■o 

14- 

i2- 

10^ 

8^ 
6^ 
4- 
2 



■ Data 
DIPHOX CTEQ6M 

■ RESBOS CTEQ6M _ 
PYTHIA YY+Yi 
PYTHIA rt 




-1 



-0.5 



0.5 



£. 0.8 
^ 0.6 

'i!; OA 

f 0.2 
S 
-0.2 
-0.4 

0. 1 

g. 0.8 
2 0.6 
a 0.4 
S 0.2 

€ 

-0.2 
-0.4 

1 

0.8 
0.6 
0.4 
0.2 
S 
-§-0.2 
"'-0.4 
-0.6 




■ Data 
PDF uncertainty 
Scale uncertainty 





-0.5 





Y.. 



0.5 





CDF II Diphoton 9.5 fb 

_ET>15,17GeV, hl<1-0, 
A R>0.4, lso<2 GeV 



Data 

NNLO 
. MCFM 
■ SHERPA 




O 

=! 0.8 
i 0.6 
5-0.4 
=! 0.2 
i -0 
n-0.2 
■g-0.4 
~-0.6 
S 1.5 



E 0.5 



0.5 
0.8 
0.6 
0.4 
0.2 
-0 
-0.2 
-0.4 



■ Data 
PDF uncertainty 
Scale uncertainly 




-0.5 





Y.. 



0.5 



Figure 10: The measured differential cross sections for ^(77), the rapidity of the diphoton system, compared with 
six theoretical predictions discussed in the text. The left windows show the absolute comparisons and the right 
windows show the fractional deviations of the data from the theoretical predictions. Note that the vertical axis 
scales differ between fractional deviation plots. 



11 



y(77) for Pt{ii) > ^(77) 



£1 
>- 



.1 

CDF II DIphoton 9.5 fb 
"ET>15,17GeV, hl<1-0, 
A R>0.4, lso<2 GeV 



■D 

e 
■a 



1.4 



1.2 



■ Data 

DIPHOX CTEQ6M 

■ RESBOS CTEQ6M 
PYTHIA YY+Yi 
PYTHIA ri 




0.5 1 
Y„forPT>M 



a. 



X 
0. 
Q 

>< 
O. 
Q 



■ Data 

■ PDF uncertainty 
Scale uncertainty 



-1 -0.5 0.5 1 

Y forPT>M 





CDF II DIphoton 9.5 fb 

"ET>15,17GeV, hl<1-0, 
A R>0.4, lso<2 GeV 



Data 

NNLO 
. MCFM 
■ SHERPA 




-0.5 



0.5 1 
Y^forPT>M 



2 

1 1.5 

^ 0.5 
? 
15-0.5 



u. 5 

S 4 

t 3 

u. 2 

^ 
E 1 

^ 

S -1 

< 

1 ^ 

°- 1 

<b 



■ Data 
PDF uncertainty 
Scale uncertainty 



r | TiTi i lT '-'-..i = i= 4^::ri.:: i f . T . ■ X,4^ i. 7r .i.| 



-1 -0.5 0.5 1 

Y forPT>M 



Figure 11: The measured differential cross sections for 1^(77), the rapidity of the diphoton system, when Pt(77) > 
M(77), compared with six theoretical predictions discussed in the text. The left windows show the absolute 
comparisons and the right windows show the fractional deviations of the data from the theoretical predictions. 
Note that the vertical axis scales differ between fractional deviation plots. 



12 



y(77) for Pt{ii) < ^(77) 



£1 
Q. 



>- 
■D 

e 
■o 



CDF II Diphoton 9.5 fb"' 
- ET>15,17GeV, hl<1-0, 
16-AR>0.4, lso<2 GeV 



14 
12 



■ Data 

DIPHOX CTEQ6M 

■ RESBOS CTEQ6M 

PYTHIA YY+Yi 
PYTHIA YY 




0.5 1 
Y„forPT<M 



a. 



1 

0.8 
F 0.6 
a. 0.4 
h 0.2 

5 
-0.2 
-0.4 

X 1 
a 0.8 

6 0.6 
S 0.4 
lb 0.2 
I 
~-0.2 

-0.4 

go.8 
= 0.6 
m 0.4 
S 0.2 
<!i -0 
1-0.2 
^-0.4 
-0.6 




■ Data 
PDF uncertainty 
Scale uncertainty 



T 




-0.5 



0.5 



for P,<M 



M- 

CDF II Diphoton 9.5 fb 

ET>15,17GeV, hl<1-0, 
-A R>0.4, lso<2 GeV 



Data 

NNLO 
. MCFM 
■ SHERPA 




-0.5 



0.5 1 
Y^forPT<M 




Y forPT<M 



Figure 12: The measured differential cross sections for 1^(77), the rapidity of the diphoton system, when Pt(77) < 
M(77), compared with six theoretical predictions discussed in the text. The left windows show the absolute 
comparisons and the right windows show the fractional deviations of the data from the theoretical predictions. 
Note that the vertical axis scales differ between fractional deviation plots. 



13 



.72 / 

T 



0.10 -CDF II Diphoton 9.5 fb"' 
^ - ET>15,17GeV, hl<1-0, 
■g a - A R>0.4, lso<2 GeV 

■a 



Data 

DIPHOX CTEQ6M 

RESBOS CTEQ6M 
PYTHIA YY+Yi 
PYTHIA YY 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Z = Et2/Eti 



1.5 
1 

0.5 

0.5 

1 

0.5 

0.5 



m J -; 

tn '■■ 

cc J.. 

m q!( 

DC 0.; 







— •— Data 

PDF uncertainty 

Scale uncertainty 


















~*~-»- ~*~ - 









































0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
z = Et2/Eti 



|10^ 

N 

lio^ 

■D 

10^ 
10== 
10 

1 

10-^ 
10-^ 



10" 



"CDF II Diphoton 9.5 fb 
ET>15,17GeV, hl<1-0, 
A R>0.4, lso<2 GeV 



Data 
NNLO 
MCFM 
■ SHERPA 




2 1 

z 0.8 

5 0.6 

0.4 

3 0.2 

1 
?-0.2 
JS-0.4 
fl-0.6 
--0.8 



1^ 0.5 
u. 

o 
2-0.5 

(0 

S -1 
S 0.8 

1 8:4 

■S-0.8 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

z = Etj/Eti 



^ — •— Data 

E- PDF uncertainty 

^ Scale uncertainty 






:irs 























0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 13: The measured differential cross sections for E'^ / E'^ compared with six theoretical predictions discussed 
in the text. The left windows show the absolute comparisons and the right windows show the fractional deviations 
of the data from the theoretical predictions. Note that the vertical axis scales differ between fractional deviation 
pfots. 



14 



Ef/E^ for Pt{ii) > M{-fj) 



Q. „ CDF II Diphoton 9.5 fb 
^10 ^ET>15,17GeV, hl<1-0, 
■g E A R>0.4, lso<2 GeV 

■a 

10- 



■ Data 

DIPHOX CTEQ6M 

■ RESBOS CTEQ6M 

PYTHIA YY+Yi 
PYTHIA YY 



10"^ 



10"' 



10" 









^^^^ 




ft 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

z = Eyj/Eyi for Pj>M 



X 

a. 
a 

>< 

0. 

a 



m 

CO 

cc 
m 

DC 



■ Data 
PDF uncertainty 
Scale uncertainty 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
z = E^j/Eti for Pt>M 



^10 

N " 

■a 

D 
■D 

10 



CDF II Diphoton 9.5 fb 

ET>15,17GeV, hl<1-0, 
E A R>0.4, lso<2 GeV 



10"^ 



10"' 



10" 



> Data 

NNLO 

— MCFM 
-■■ SHERPA 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

z = Eyj/Eji for Pj>M 



< 
a. 

cc 

CO 

< 

0. 

cc 

CO 



7 
6 
5 
4 
3 
2 
1 


■i 

6 
5 
4 
3 
2 
1 


■n 

61 

5 

4 

3 

2 

1 



-1 



■ Data 
PDF uncertainty 
Scale uncertainly 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 14: The measured differential cross sections for when Pt(77) > Af(77), compared with six 

theoretical predictions discussed in the text. The left windows show the absolute comparisons and the right 
windows show the fractional deviations of the data from the theoretical predictions. Note that the vertical axis 
scales differ between fractional deviation plots. 
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Figure 15: The measured differential cross sections for when Pt (77) < M{'jj), compared with six 

theoretical predictions discussed in the text. The left windows show the absolute comparisons and the right 
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Figure 16: The measured differential cross sections for the cosine of the polar angle of the leading photon in 
the CoUins-Soper frame compared with six theoretical predictions discussed in the text. The left windows show 
the absolute comparisons and the right windows show the fractional deviations of the data from the theoretical 
predictions. Note that the vertical axis scales differ between fractional deviation plots. 
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Figure 17: The measured differential cross sections for the cosine of the polar angle of the leading photon in the 
CoUins-Soper frame, when Pt{ii) > M^jj), compared with six theoretical predictions discussed in the text. The 
left windows show the absolute comparisons and the right windows show the fractional deviations of the data from 
the theoretical predictions. Note that the vertical axis scales differ between fractional deviation plots. 
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Figure 18: The measured differential cross sections for the cosine of the polar angle of the leading photon in the 
CoUins-Soper frame, when Pt{ii) < -^^(77)1 compared with six theoretical predictions discussed in the text. The 
left windows show the absolute comparisons and the right windows show the fractional deviations of the data from 
the theoretical predictions. Note that the vertical axis scales differ between fractional deviation plots. 



19 



m - m 



S 12 

< 

■D 



D 

■a 



n \ \ I ^ \ I \ \ I \ \ r 
CDF II DIphoton 9.5 fb'' 
ET>15,17GeV, hl<1-0, 
A R>0.4, lso<2 GeV 



10 



-3 



I \ \ I \ ^ I \ r 

■ Data 

DIPHOX CTEQ6M 

■ RESBOS CTEQ6M 

PYTHIA YY+Yi 
PYTHIA YY 



-fj^ ..... '"K't. 



0i ■:■ 



-2 



-1 



At) 



1 

0.8 
0.6 
0.4 
0.2 

■0.2 

0-1 
0.8 
0.6 
0.4 
0.2 

■0.2 

°-1 
0.8 
0.6 
0.4 
0.2 

■0.2 
■0.4 




■ Data 
PDF uncertainly 
Scale uncertainty 




-2-10 12 3 

At| 



"1 — \ — \ — I — ^ — \ — I — \ — \ — I — \ — \ — r 
CDF II DIphoton 9.5 fb'^ 
ET>15,17GeV, ^l|<1.0, 
A R>0.4, lso<2 GeV 



■ Data 

PYTHIA YY+Yi 

MCFM 
■ SHERPA 




1 

0.8 
0.6 
0.4 
0.2 

■0.2 
0.4 

1.2 
1 

0.8 
0.6 
0.4 
0.2 


■0.2 
0.4 
0.8 
0.6 
0.4 
0.2 

-0 
-0.2 
-0.4 




<^-J:±. ^ ^ ^ 



-4- 





At] 



Figure 19: The measured differential cross sections for r]i — 772, the diffrence in photon pseudorapidities, compared 
with six theoretical predictions discussed in the text. The left windows show the absolute comparisons and the 
right windows show the fractional deviations of the data from the theoretical predictions. Note that the vertical 
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Figure 20: The measured differential cross sections for logio (^^(77)/-^" (77)) compared with six theoretical pre- 
dictions discussed in the text. The left windows show the absolute comparisons and the right windows show the 
fractional deviations of the data from the theoretical predictions. Note that the vertical axis scales differ between 
fractional deviation plots. 
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Figure 21: The measured differential cross sections for Ai?(77) where Ai?(77)^ = A<I>(77)^ + {rji — 772)^ compared 
with six theoretical predictions discussed in the text. The left windows show the absolute comparisons and the 
right windows show the fractional deviations of the data from the theoretical predictions. Note that the vertical 
axis scales differ between fractional deviation plots. 
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Figure 22: The measured differential cross sections for the photon Et (two entries per event) compared with 
six theoretical predictions discussed in the text. The left windows show the absolute comparisons and the right 
windows show the fractional deviations of the data from the theoretical predictions. Note that the vertical axis 
scales differ between fractional deviation plots. 
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Figure 23: The measured differential cross sections for the photon rj (two entries per event) compared with six 
theoretical predictions discussed in the text. The left windows show the absolute comparisons and the right windows 
show the fractional deviations of the data from the theoretical predictions. Note that the vertical axis scales differ 
between fractional deviation plots. 
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